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Abstract 

The quark mass dependence of hadrons is an important input for lattice calculations. We 
investigate the light quark mass dependence of the binding energy of the Y(3872) in a finite 
box to next-to-leading order in an effective field theory for the Jl( 3872) with perturbative pions 
(XEFT). At this order, the quark mass dependence is determined by a quark mass-dependent 
contact interaction in addition to the one-pion exchange. While there is only a moderate sensitivity 
to the light quark masses in the region up to twice their physical value, the finite volume effects 
are significant already at box length as large as 20 fm. 

PACS numbers: 14.40.Pq, 13.75.Lb, ll.30.Rd 
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I. INTRODUCTION 


The discovery of the X(3872) in 2003 by the Belle Collaboration [1] with its conhrma- 
tion by the CDF collaboration shortly after [2], was the hrst of a series of discoveries of 
charmonium-like hadrons [3]. Its decays into [4] and J/'ilJ7r~^7r~7r^ [5] with a ra¬ 

tio of branching fractions close to one, indicate a large isospin violation in X decays and 
make an interpretation as a conventional cc state implausible. The description as a loosely 
bound D^D*^ S'-wave hadronic molecule with even charge parity [6-11], on the other hand, 
naturally explains the proximity of the X mass to the D^D*^ threshold and the quantum 
numbers = I’*"’'' [12, 13]. 

In the molecular picture, the binding energy. Ex, is determined by the masses of the 
X(3872), the and the D*^ meson, Mx, mo and mo*, respectively. Using the latest 
values from the review of particle properties [14], it reads 

Ex = rriD* + rriD — Mx = (0.11 ± 0.21) MeV. (1) 

A natural energy scale is set by the one-pion exchange, m‘^/{2MDD*) ~ 10 MeV, where 
Mr,]:,* is the reduced mass of the and D*^ mesons and the neutral pion mass. The 
binding energy is small compared to this natural energy scale and hence, the X(3872) 
displays universal properties. 

Braaten and Kusunoki exploited this universality in a series of papers on the X using 
effective held theory methods [10]. They obtained various predictions for production am¬ 
plitudes, decays, formation, and line shapes of the X(3872) (see Ref. [15] for a review). In 
Ref. [16], the binding energy of the X(3872) was calculated in a pionless effective held theory 
using constraints from heavy-quark symmetry. The inhuence of three-body DDtt interac¬ 
tions on the properties of the X(3872) was found to be moderate in a Faddeev approach 
[17]. Finally, we note that universality also determines the interactions of the X(3872) with 
neutral D and D* mesons [18]. The corrections to universality can be calculated systemati¬ 
cally using an ehective held theory for the X with explicit pions, called XEFT, which was 
developed by Fleming, Kusunoki, Mehen and van Kolck in 2007 [19]. They applied XEFT 
to calculate the partial decay width F [X —)■ D°D°7r°] at next-to-leading order (NLO) in the 
XEFT power counting. Later, their work was extended to describe hadronic decays of the 
X(3872) to the Xcj [20]. In Ref. [21], it was pointed out that XEFT can also be extended 
to systems with an additional pion with energies close to the D*D* threshold. Finally, a 
Galilean-invariant formulation of XEFT was introduced in Ref. [22] to exploit the fact that 
mass is very nearly conserved in the transition D*^ —)■ 

Whereas the X(3872) is appealing for ehective held theory approaches particularly be¬ 
cause of its unnatural size, its large extent poses a severe technical problem for lattice 
simulations. However, a quenched lattice calculation supported the = I’*"’'' hypothesis 
[23] before the LHCb experiment hnally settled the X quantum numbers [13]. A full lattice 
QCD study was hrst performed by Prelovsek and Lescovec in 2013 [24]. In this calculation a 
candidate for the X(3872) about (11 ± 7) MeV below the D^D*^ threshold has been found. 
The authors used light quark masses at about four times the physical value and a spatial 
lattice size of only L = 2 fm. The typical length scale of the X can be estimated from the 
jjOD*^ R-wave scattering length at leading order: = 1 /i/2Md^TEx [25]. For the X, this 

implies a* > 5 fm ^ L, such that large chiral and hnite volume ehects are expected for the 
calculation of [24]. Two recent lattice studies [26, 27] use similar volumes and pion masses 
such that these problems persist. While the quark mass dependence has been addressed in 
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Refs. [28-31], a calculation of the finite volume corrections to observables of the X(3872) is 
still outstanding. 

For two particles in a finite volume, Liischer has developed a framework to determine 
bound-state and scattering observables from hnite volume energy levels [32, 33]. An equiv¬ 
alent pionless effective held theory approach for two nucleons on a lattice using the power 
divergence subtraction (PDS) scheme was presented in [34]. In the quark mass range where 
the D*^ is unstable, however, the D^D*^ system has on-shell intermediate states 

and signihcant three-body effects are expected. In the last years, several attempts have 
been made to obtain a better understanding of the three-body system in a box. The modi- 
hcation of three-body bound states in a cubic volume was investigated in pionless effective 
held theory [35, 36]. In Refs. [37, 38] and [39], it was shown that hnite volume observables 
are determined by the inhnite volume S'-matrix elements as it is the case for the two-body 
system in a box. There are also some explicit calculations available for systems possessing 
three-body intermediate states. The ai(1260) resonance in a cubic box, generated by np 
scattering and accounting for the p meson’s self energy, has been addressed in [40]. Further, 
in Ref. [41], topological ehects for bound states in a moving frame have been considered. 
An analytical expression for the hnite volume energy shift of three identical bosons in the 
unitary limit was derived in Ref. [42]. 

In this work, we investigate hnite volume corrections to the binding energy of the X (3872). 
We present explicit expressions to extrapolate results for the binding energy of the X(3872) 
from the hnite to the inhnite volume and from unphysical to physical quark masses. We 
consider both, quark masses for which the D*^ is stable and quark masses for which it can 
decay into . Furthermore, we use the quark mass dependence synonymous to the pion 
mass dependence because of the Gell-Mann-Oakes-Renner relation [43]: 

ml =-{mu +rnd){^\uu +dd\^)/P , (2) 

where / ~ 130 MeV is the pion decay constant, and are the light quark masses, and 
(0|-un|0) = (0|(i(i|0) = (—283(2) MeV)^ is the light quark condensate in the MS scheme at 
2 GeV [44]. 

The paper is organized as follows: In Sec. II, we briehy review XEFT in the inhnite 
volume. A strategy how to obtain shifts for the binding energy due to higher order con¬ 
tributions is presented in Sec. III. The hnite volume amplitudes are calculated in Sec. IV 
and explicit expressions for the binding energy in dependence on the box size are given. In 
Sec.V, we discuss our results for the chiral and hnite volume extrapolations of the X(3872) 
binding energy. Finally, we summarize our conclusions and present an outlook on future 
work in Sec. VI. 


II. REVIEW OF XEFT IN THE INFINITE VOLUME 

We briehy review the results for the binding energy in the inhnite volume, following 
our analysis in [30]. The underlying ehective theory we are using in order to describe the 
X, called XEFT, was derived by Fleming et. al. in [19] starting from heavy-meson chiral 
perturbation theory. In XEFT, regarding the X(3872) as a D^D*^ A-wave hadronic molecule, 
£) 0 , jj*o .^0 ggifig treated non-relativistically. Gharged 77^) mesons are 

integrated out and do not contribute to the order we are working at. Moreover, the 
and J/'^ca channels can also be integrated out [19]. Utilizing XEFT, we can evaluate the 
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center-of-momentum D^D*^ S'-wave scattering diagrams and eventually extract the binding 
energy of the X. 

The XEFT Lagrangian reads 


— I t I I I 

\ irriD* I \ zrriD I 

+D'> (ido + ^ D + (ido + ^ D + (id^ + 

\ ItTId* I \ IlTljj I \ 2171^ 

' ^ ^ (^DD^ ■ + DW ■ + h.c. 


^/2f \/2m7r ' 

y {DD + DDy ■ {DD + DD) 

%{dd + dd)'-[ 


+^(DD + DDy- 

(dd + dD) 


^ + dV'^D^ + h.c. 

(DD + DD) +■■■ . 


( 3 ) 


Here, v = (ht v —tnv)/(tn + r^) is the Galilean invariant derivative, In (nt) the mass of 
the left- (right-) hand held and the ellipsis denote higher-order interactions. The masses for 
the pion, the D^- and Zi)*°-meson are labeled by ttItt, niD and mn*, respectively. Furthermore, 
g is the D meson axial coupling constant, / the pion decay constant and A = tud* — mo the 
hyperhne splitting of the D mesons. The mass scales fi and <5 are dehned as 
and S = A —m^. The coupling constants Cq, C 2 and D 2 are discussed below. Note that this 
Lagrangian has no exact Galilean invariance. The D*^D^7r^ coupling is given by the leading 
term in the chiral limit and only the leading terms in expansion in m.,^/m£) and (5/m^ are 
kept in the calculated observables. 

In its structure, XFFT is similar to the Kaplan-Savage-Wise theory (KSW) for nucleon- 
nucleon (NN) scattering, which makes use of the power divergence subtraction scheme 
(PDS) [45]. PDS has proven to be well suited for systems with an unnaturally large scattering 
length. In the KSW counting, the pion exchanges are included perturbatively. Although 
the perturbative treatment of the pion exchanges was shown to fail in the NN sector at 
NNLO because of large contributions from the nuclear tensor force [46], it is expected to 
work well for D^D*^ scattering in XEFT due to a signihcantly smaller expansion parameter 
[19]. In XEFT, the mass scale /x, the momenta of the mesons and the pions as well 
as the binding momentum are all counted as order Q, which dehnes the typical momentum 
scale of XEFT. At leading order (LO), that is Q~^ in XEFT power counting, there is one 
contact interaction with coupling constant Cq. Loop integrations contribute a factor 
and propagators a factor Q~^. Therefore a loop {Q^) consisting of two propagators (each 
Q~^) and a LO contact interaction (Q“^) is counted as order . Appending such a loop 
to any diagram leaves the order of the diagram unchanged. On the one hand, this implies 
that the LO contact interaction has to be resummed to all orders. On the other hand, the 
higher order contributions have to be dressed in all possible ways by LO amplitudes. The 
NLO amplitude is of order and three more interactions have to be considered: two NLO 
contact interactions with coupling constants C 2 and D 2 and the D*^D^7r^ coupling. 

Further, XEFT explicitly accounts for the hnite decay width of the D*^ meson. Whereas 
nucleons, regarding strong interactions only, are stable, the D*^ can decay into D^ti as 
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long as the hyperfine splitting of the D*^ and the meson is greater than the pion mass. 
Thus additional diagrams have to be included which lead to the emergence of infrared (IR) 
divergences. In order to cure this pathological behavior, we resum pion contributions to the 
propagator of the D*^ meson [30]. Inserting these resummed instead of bare D*^ propagators 
in all D^D*^ loops, we obtain IR hnite amplitudes. For more details on infrared divergences, 
the XEFT power counting and the derivation of the Lagrangian, we refer to [19] and [30]. 

In the Galilean-invariant version of XEFT from [22], the width from decays of the D*^ into 
both and D^'-y is included and the Galilean invariance is exact. It strongly constrains 

the form of the ultraviolet divergences in the theory such that no expansion in m.,^/m£) 
is required. A covariant formulation of a non-relativistic effective held theory describing 
Goldstone boson dynamics can be found in Ref. [47]. 




VOS 


— fS + 


+ X 


FIG. 1: Self-energy graph with pion contribution and counterterm for the D*^. Single, double and 
dashed lines represent , D*^ and pion propagators, respectively. The cross on the right hand 
side indicates the insertion of the counterterm. 

We begin with the derivation of the resummed D*^ propagator, utilized in all calculations 
for the D^D*^ system. First, we consider the D*^ self-energy shown in Fig. 1. For the 
counterterm, we use the on-shell renormalization scheme. It ensures that the real part 
of the propagator’s pole position is at the on-shell point po = P^/2m£)», with po being the 
energy and p the momentum of the D*^. Using PDS^, the bare self-energy reads 

^ 2^ 

where A is the PDS renormalization scale. The counterterm is chosen such that it cancels 
the second term in parentheses, which is real valued, analytic in the quark masses and 
proportional to the PDS renormalization scale A. The hrst term is purely imaginary as long 
as the D*^ can decay into i.e. for pion masses smaller than the hyperhne splitting. 

In this case it induces a hnite decay width. For pion masses greater than the hyperhne 
splitting, it is real valued and the self-energy implies a hnite mass shift for the denoted 
by . In summary we have 


= 


Am^,* 


W 


24 : 71 / 




0 , 


247r/' 


2X^ 


m. 


< A, 
> A. 


( 5 ) 

( 6 ) 


We point out that the mass scale p is purely imaginary above the decay threshold and 
hence A^^, is real valued for all pion masses. The D*^ propagator can now be calculated 


^ Note that the additional term occurring in PDS, proportional to the renormalization scale A, is subtracted 
again due to the use of on-shell renormalization scheme. 
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according to Fig. 2 and we obtain after resummation 
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Po ~ p‘^/2m£)* + SO® + ie 


( 7 ) 
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FIG. 2: Fully resummed D*^ propagator with the insertion of the self-energy Bare and full 

propagator are represented by double and thick lines, respectively. 


We proceed with the evaluation of the D^D*^ scattering diagrams. Since we consider 
S-wsive scattering, the total angular momentum of the D^D*^ system is determined by the 
D*^ meson’s spin. Denoting the polarization vectors of the incoming and outgoing D*^ by 
Si and £*, respectively, it turns out that all amplitudes factorize and we can write 

( 8 ) 


using spin indices i and j. For the calculation of the binding energy, it is sufficient to take 
a look at the scalar amplitudes A. At leading order, there is one contact interaction with 
coupling constant Gq. According to XEFT power counting, it has to be resummed to all 
orders, as depicted in Fig. 3. Using PDS to renormalize the linear divergence of the loop 
integral in Fig. 3, 


h = 


1 

(27r)^ |k|^ -|- T]'^ 


PDS 1 /. X 


(9) 


with the energy-dependent quantity r] = 2 Moo*(U -|- S^®) — ie, we obtain 

2TTi 1 27ri 1 


iA-i = 


M, 


DD* 


2n 


The LO coupling constant occurs in the definition of 7 

271 


dvr/o (p 2 ) Mdd* -7 + ^' 


7 = 


Mdd*Co{A.) 


A. 


( 10 ) 


( 11 ) 


The LO amplitude has a pole at —E = = 7 ^/ 2 M£)£)* -|- and we identify the binding 

energy with the real part of . Thus, 7 is the LO binding momentum, and Cq in Eq. 
(11) has to cancel the dependence on the PDS renormalization scale A and fix the binding 
energy for physical quark masses to the experimentally measured value. 

At NLO, there are three more interactions, which are included perturbatively. Two 
contact interactions, with coupling constants C 2 and D 2 , and the one-pion exchange (OPE). 
Following XEFT power counting, the LO amplitude has to be appended in all possible ways 
to the NLO interactions. We end up with the scattering diagrams shown in Fig. 4. The 
results for the S'-wave projected scalar amplitudes are given in Appendix A. 

After the inclusion of NLO contributions, the renormalization condition, when fixing the 
binding energy at its experimentally measured value for physical quark masses, implies a 
relation between the LO and the NLO coupling constants. 
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FIG. 3: Leading-order diagram for D^D*^ scattering. The single and double lines represent 
and D*^ mesons, respectively. The thick, internal line indicates a fully resummed D*^ propagator 
depicted in Fig. 2. 
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FIG. 4: Next-to-leading-order diagrams for D^D*^ scattering. We use the same notation as in 
Figs. 1, 2 and 3. 

Since C 2 and D 2 are unknown, we estimate natural ranges and vary the coupling constants 
to determine the error band of the binding energy. We rewrite the coupling constants as 

C 2 = = Cj(C„)", (12a) 

= (C'ot, (12b) 

where i? is a renormalization constant given in Appendix A and the superscript “ph” indi¬ 
cates that a quantity, here the mass scale /i, is evaluated at physical quark masses. Further¬ 
more, compared with the pure contact theories, tq can be identihed with the effective range 
in the pionless theory. We follow [19] and [30] and use their estimates tq E [0,1/lOOMeV] 


7 
























and d 2 G [—1,1]. In the future, it should be possible to determine C 2 and D 2 from lattice 
calculations. 


III. STRATEGY FOR EXTRACTING THE BINDING ENERGY TO NLO 


For an unstable D*^ meson, the OPE potential is oscillatory and not Yukawa-like [48]. The 
effective range expansion breaks down at NLO and the effective range is not dehned. Hence, 
the binding energy can not be extracted from effective range parameters. In this section, we 
present an alternative method to access the binding energy up to NLO, employing the two- 
body scattering amplitudes, regardless of whether examining stable or unstable particles in 
the hnite or inhnite volume. First, we note that the sum of the NLO scattering amplitudes, 
Aq \ • ■ ■ , can be collected in powers of the LO amplitude 


Ao — Aq ^ Aq^^ — -So -f siA_i -|- S2 (A_i)^. 


Furthermore, we expand the LO amplitude around the LO pole position 


A_i = 


^-1 

E + HLO 


+ ••• 


1 


(13) 


(14) 


where the dots denote terms being hnite at E = and Z_i is the residue 




i 


d 1 

^iA_i 


_B=-BL0 


— {MDD*y 1 
2n 7 


(15) 


Accordingly, the full amplitude up to NLO, expanded around the LO pole position, can be 
written as 


A 


A_i -l- Ao — 


Z_i + SiZ_i S2{Z_iy 

E + B^o ^ (E + + ■ ■ ■ • 


(16) 


Moreover, we consider a generic, non-perturbative expression for the amplitude with shifted 
pole position, B = B^^ + AB, and shifted residue, Z = Z_i -|- AZ, and expand it around 
the LO pole position 


Z , _ Z_i + AZ ZAB 

E + B ^ E + B^° ~ (E + B“)2 ^ ' 


(17) 


Utilizing expressions (16) and (17), the NLO shifts for the residue, Z_i, and the LO pole 
position, B^^, can be read off by equating coefficients 


^^NLO ^ 

AHNLo = ^2 


1 + <§1 


Z-i -)■ —S 2 Z- 1 , 


(18a) 

(18b) 


where we already used a partial NNLO cancellation in Eq. (18b), which is described in 
Appendix B in more detail. For a pure contact theory, a comparison to an approach, where 
the NLO coupling constants are resummed to all orders, can be found in Sec. V. In the 
following, we apply this strategy to extract the binding energy of the X(3872) in a hnite 
volume. 



IV. FINITE VOLUME CORRECTIONS TO THE BINDING ENERGY 


We consider the D^D*^ system in a box with side length L and periodic boundary condi¬ 
tions. The allowed lattice momenta are then given by integer vectors times 2ti/L. Integrals, 
occurring in calculations for the binding energy in the inhnite volume, have to be replaced 
by discrete sums over the quantized lattice momenta. Since we are interested in D^D*^ 
bound states with even parity, we expect that the binding energy acquires a positive shift. 
We distinguish between two different regions. One, where the D*^ is unstable, i.e. for pion 
masses < A and a second region where the D*^ is stable, i.e. > A. Whereas in the 
hrst case explicit XEFT calculations have to be carried out due to three-body DDn inter¬ 
mediate states, in the latter case one can alternatively use a two-body approach introduced 
in [34], which serves as a consistency check. All quantities which differ in the hnite volume 
are tagged by a superscript L. 


A. The LO amplitude 


Let us begin with the explicit XEFT calculations, which can be utilized in both regions. 
To LO, the amplitude in the hnite volume reads 


iA^i 


2Tii 


^DD* 


27r 




Co 


’ 


(19) 


where the hnite volume quantity Iq is given by 





n G 


( 20 ) 


Since the mass of the D*^ meson obtains a shift in the hnite volume as given below, the 
reduced mass is dependent on the box size, too. Like the loop integral (9) in the inhnite 
volume, Iq is linearly ultraviolet divergent since short-distance properties of the theory 
remain unchanged in the hnite volume. We regularize following [34]: hrst we introduce a 
sharp momentum cut-oh. A, for the sum and then add and subtract the inhnite volume loop 
integrals evaluated at zero energy. One of the loop integrals is regularized using PDS, the 
other one using a momentum cut-oh, which coincides with the cut-oh in the sum. Finally, 
the limit A —)■ cx) is taken. We obtain 


PDS 


» lim 

A—^CXD 


. |k|<A 

— V — 

1 

J 

f d^k e{x 

- 

■|k|) 

L3 ^ |k| 

J 

{2nf 


2 



( 21 ) 


with E = Plugging (21) into (19) and using the dehnition for the LO binding 

momentum in the inhnite volume (11), we acquire 


iA^i = 



( 22 ) 
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with 


S{x) 


lim 

—^oo 


|n|<A„ ^ 

772 47 rA„ , 


(23) 


where A„ = LA/(27r). The energy levels of the D^D*^ system to LO in the hnite volume 
can be determined from Eq. (22). Note that they are fully determined by the inhnite 
volume quantity 7 .^ Here, we are interested in the solution with negative energy, i.e. the 
solution which approaches the inhnite volume LO binding energy for L ^ 00 . We denote 
the corresponding LO binding momentum by 7 '^, dehned by 



(24) 


B. The D*^ self-energy and mass shift 


We proceed with the D*^ self-energy and mass shift. The calculation is carried out 
similarly to the one of the LO amplitude. Using PDS and cut-off regularization we obtain 
for the bare self-energy 




2Mf 


L\2 




L\2 





(25) 


Independent of the pion mass, we do not receive any imaginary contributions for the bare 
self-energy in the box. However, since the hnite volume itself cuts oh low frequency modes, 
we do not expect the occurrence of any infrared divergences. 

Again, we use the on-shell renormalization scheme and subtract the second term propor¬ 
tional to the PDS renormalization scale A. Hence, the counterterms in the hnite and inhnite 
volume coincide up to corrections to / and /x. The shift for the D*^ mass is diherent though. 




2Mf 


L\2 




2n 


(26) 


Note that even for physical pion mass, the D*^ meson in a box receives a hnite mass shift. 


C. NLO corrections to the binding energy 

Now, we implement the corrections due to the NLO amplitudes. In analogy to the inhnite 
volume we hud for the NLO contact interactions, i.e. the amplitudes and 

*“^ 0 ( 1 ) = (27a) 

WoJ 

i.4S,v, = (27b) 


^ We expect that even if three particle intermediate states exist that is after the inclusion of NLO contri¬ 
butions and for < A, finite volume observables are still determined by the infinite volume S'-matrix 
as demonstrated in [37-39]. 
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For the pion exchange diagrams, we do not project onto the S'-waves. Whereas the inhnite 
volnme is rotationally invariant, the lattice is only invariant nnder transformations of the 
cnbic gronp. In principle, it is possible to decompose qnantities transforming according to 
an irredncible representation of the cnbic gronp into spherical harmonics [49, 50]. However, 
we keep the sums over integer vectors, since convergence of the partial wave expansion is 
not certain. The OPE amplitude ^o(ii) is then given as 


'i A ^ 




2{f^y 

with the incoming (outgoing) relative momentum £ {£'). For ^o(iii) we hnd 




L 


ig 


271 2(/^)2 

where the quantity is dehned as 


(f. (¥)')+p ^ ^ 


(28) 


(29) 


c-(ni) 


(m, x) 


lim 

Xn —^OO 


|n|<A„ 

E 


fn + m) £* ■ (n + m) 


n 


m 


n + m 


X 




—47rA^ 
3 


(30) 


The amplitudes and ^o{iii)ij imply a coupling between channels with different angu¬ 

lar momentum. Considering the Ai representation of the cubic group, the lowest angular 
momenta coupled are with I = 0, 4, 6, 8, • • ■. On the other hand, for the amplitude ^o(iv)ij 

we can use a tensor decomposition and it appears that vAg^jy^^^. = ^p-Ag^jy). A detailed 
derivation is given in Appendix C. We obtain for the scalar amplitude 





lim 

An—^OO 


|n|,|n^|<An 

E 

r» n' 




I n -|- n' I 


27r" 



(31) 



Due to the coupling between different angular momenta, the amplitudes and 

in the hnite volume do not factorize into a scalar amplitude and a function of the incoming 
and outgoing D*^ mesons’ spins, in particular A This implies a non¬ 

trivial dependence of the coefficients sq and si in Eq. (13) on the polarization vectors e* 
and £j and hence of the NLO shift for the held strength renormalization constant, 

However, since the amplitudes and do factorize^ and therefore S 2 , it 


^ The factorization takes place since the amplitudes Agj-jY^y and Agj-y^y contain the momentum 

and hence angular independent LO amplitude A-i on both sides and thus are momentum and angular 
independent by themselves. 
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is sufficient to consider the scalar amplitudes ^o(i)’ “^^(iv) calculate the shift 

for the binding energy. 

The dependence of the loop integrals on the PDS renormalization scale A is the same as for 
the inhnite volume and accordingly the NLO coupling constants coincide with the ones given 
in Eqs. (12a) and (12b) up to hnite volume corrections to scales large compared to Q like 
for example mn* or /. This corresponds to a multiplicative renormalization scheme where 
loop integrals are regularized separately. Again, the error bands are obtained by varying 
the coupling constants within their natural ranges. For the binding energy we employ the 
results of the previous section. The quantities Z_i and S 2 have to be reevaluated in the box. 
We hnd for the residue 


{z\r^ = 


d 




\ 27r 


^=-■5^0 


[ dE 

where B^q = ( 7 -^)^/ 2 M^^* and 

S"(a;) = dxS{x) 


271 


tiL 


S' 


E 


|2 \ 
n| — X) 


2 • 



(33) 


(34) 


For the coefficient we obtain, already inserting the redehnitions of the coupling constants. 


4 = C2(74^ + 


6 (/ 


L\2 


X 


L\2 


(7 - A)^ + (4) 


271 

- 6.2 + log 


A 


L,ph 


14 


1 

+ 2 



where we used Fq. (24) for the first term in parentheses. 


D. Validity range of XEFT in the box 


In the inhnite volume, the range of applicability of XFFT is constrained by two demands. 
On the one hand, we require that pions can be included perturbatively, determining the 
boundary for large pion masses. On the other hand, treating pions non-relativistically 
settles the low boundary. In summary, we have in the inhnite volume 0.98(m4)^ ^ 
m2 <2(m4)2 [30]. 

However, for three particles in the hnite volume, singularities occur as soon as three-body 
propagators can go on-shell, a behavior which has already been investigated e.g. in [37] and 
[38]. In XFFT, this manifests in the last term of Fq. (35). For pion masses smaller than 
the hyperhne splitting, where the D*^ —)■ D^tt decay channel is open, possesses 

singularities for values of (L/i'^/ 27 r )2 being the absolute value of an integer vector squared, 
greater or equal than one. Since the D*^ —)■ D^ti decay proceeds via a P-wave interaction, 
^2 is hnite for (L/i^/27r)2 = 0. So for certain values of m^r and L, the perturbative treatment 
clearly fails. To obtain a region of validity for XEFT in dependence on the volume and the 
pion mass, we take a look at the quantity 
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which explicitly accounts for the singularities of for < A and approaches the 

inhnite volume XEFT expansion parameter for > A and L ^ oo. The second term 
in Eq. (36) ensures that is finite for /i'^ —)■ 0. A density plot is shown in Fig. 5. We 



FIG. 5: Density plot for the expansion parameter White regions belong to values of L and 
where the perturbative inclusion of pions breaks down. For pion masses close to the physical value, 
there are several narrow bands of hyperbolic shape due to on-shell three-body propagators. 

restrict our analysis on regions where 6 ,^ < 0.15 such that it is small enough to compensate 
for unnaturally large NNLO coefficients of similar size as in KSW. For physical pion mass 
it follows that a perturbative treatment of pions is justified for 5 fm < L < 20 fm. We 
point out that the NLO parameters C 2 and d 2 coincide in the finite and infinite volume and, 
once determined from lattice calculations, one can utilize the infinite volume formulas to 
extrapolate to L ^ oo. 


E. Effective range expansion for large mn 

The preceding analysis is valid for all pion masses. Now, we focus on the region where 
the D* is stable, i.e. for > A. Here, we can apply the effective range expansion for 
the infinite volume amplitude, analytically continue it to negative energies and apply the 
procedure established in [34]. We introduce the S'-wave scattering phase shift Ss{p), which 
is related to the inhnite volume scattering amplitude by 

27r 

pcot6s{p) = tp+— -j (37) 

and apply the effective range expansion 

pcot Ss{p) = —- + H-. (38) 

dg L 

The quantities and are known as S'-wave scattering length and S-wave effective range, 
respectively. In the pionless theory, l/a^ = 7 and = tq. However, including pions leads 
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to corrections of NLO, which can be determined by expanding the inverse inhnite volume 
scattering amplitude in powers of p. Equating coefficients yields the following expressions 
for tts and 
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Let us briefly consider an effective theory in the inhnite volume, where pion interactions 
are not included explicitly but via modihed LO and NLO coupling constants with similar 
renormalization condition as in (12a) with tq replaced by in Eq. (39b). The criteria for a 
bound state follows from Eq. (37) and reads, applying Eq. (38) and neglecting higher order 
shape parameters 

— + - 7* = 0, (40) 

Clg / 

where 7* is the binding momentum including NLO contributions. Then the binding energy 
up to NLO is given by 
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(41) 


Using the same effective theory but now including NLO corrections using the strategy de¬ 
scribed in Eqs. (13) through (18), we obtain the same result as in the second line of Eq. (41) 
except that no terms of order r^/a^ occur. Hence, as long as the 5'-wave scattering length is 
signihcantly larger than the S'-wave effective range, both methods deliver consistent results. 

Using a pionless effective held theory in the hnite volume, the amplitude can be calculated 
in analogy to Sec. IV A and the result is given by Eq. (22) with 7 replaced by —pcot5s{p). 
The criteria for a bound state looks similar to (24) and (40) (cf. [34]) 


1 

Qjg 





0 . 


(42) 


Here yf" is the hnite volume binding momentum including NLO corrections. Eq. (42) 
approaches Eq. (40) in the limit L —)■ cx). As for the inhnite volume, we expect that the 
results from the two diherent methods agree as long as 3> r^. 


V. RESULTS 

In order to determine the hnite volume and quark mass dependence of the binding energy, 
we hrst consider the extrapolations for the pion decay constant, the D meson axial coupling 
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constant and the and D*^ meson masses, respectively. A superscript (0) denotes the 
chiral-limit value of a quantity. For the chiral extrapolation of the pion decay constant we 
use the results from [51] 
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with the low-energy constant U = 4.4 and = 124 MeV, corresponding to = 132 MeV 
[51, 52]. Further, we use the lattice results from [53] for the quark mass dependence of the 
D meson axial coupling constant 




1 + 2c/W^ 




■m^ log - am^ 

hlat. 


47r2j(o)2 

where the parameters are given as [53] 

= 0.46, a = —0.16 GeV“^, piat = 1 GeV. 


(44) 


(45) 


The D meson axial coupling constant does not receive any corrections in the hnite vol¬ 
ume. For the pion decay constant we employ the results given in [54] obtained from chiral 
perturbation theory to one loop 


^TT 1 A'ldnl m,r-h) 
27r/2 ttL |n| 


(46) 


with Ki being the modihed Bessel function of second kind. The chiral and hnite volume 
extrapolations for the and meson masses can be summarized as 


L ph I / 2 ( ph\2\ 

= rriD = + ^(m^ - [ml ) ), 


m 


ph ' 
D 




rrijy, = m 


D* 


m 


ph 

D* 


7Tl * 


(47) 

(48) 


with hi = 0.42 [55] and given in Eq. (26). 

In Fig. 6, we plot the dependence of the binding energy on the side length of the box, L, 
to compare the two approaches described in chapter IV. The pion masses are hxed at values 
of = 145 MeV and = 160 MeV, respectively. The inhnite volume results are shown 
by solid lines. The upper bound corresponds to values of the NLO parameters of ^2 = 1 
and ro = 0 and the lower bound to ^2 = and tq = 0. This results in the maximum band 
width. In the hnite volume, the parameter values for the lower bound are the same but the 
upper band belongs to d 2 = 1 and tq = 0.01/MeV in order to maximize the error band. 

Whereas the lower bounds and central values coincide well using the two diherent strate¬ 
gies and deviations are clearly smaller than the NLO shifts, there is some discrepancy for 
the upper bounds. This can be understood from the considerations in Sec. IV E. Results 
are consistent as long as the S'-wave scattering length is much larger than the S'-wave ehec- 
tive range. For ^2 = 1 and tq = 0.01/MeV however, and are similar in size and the 
error induced by approximating the root in (41) is > 10%. This is in the order of the NLO 
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FIG. 6: Comparison of the two methods to obtain the binding energy to NLO described in chapter 
III and IV. We keep the pion mass hxed at values of 171-,^ = 145 MeV (left) and rrin = 160 MeV 
(right). The binding energy in the infinite volume is represented by solid lines. The band is 
acquired by varying the NLO parameters d 2 and rg within their natural ranges. The thick, central 
curve belongs to ^2 = 0 and rg = 0. The finite volume results correspond to crosses for the explicit 
XEFT calculations and to empty squares for the results obtained from an effective range expansion. 
The central crosses and squares belong to d 2 = 0 and tq = 0. 


corrections and explains the deviation for the npper bonnds in Fig. 6. We point out that 
effective range and scattering length are of comparable magnitude for a very limited range 
of the NLO parameters only. 

So far we looked at pion masses above the hyperhne splitting of the H-mesons. We now 
consider the region where the D*^ can decay into D^tt. The binding energy in a hnite volume 
for physical pion mass is depicted in Fig. 7. We plot for box lengths between 5 and 20 fm 
where the expansion parameter in Eq. (36) is clearly smaller than 0.15 as can be read off 
from Fig. 5. 

The result for d 2 = 0 and ro = 0 is not shown as it almost coincides with the lower bound. 
This can be understood by noting that the only difference of the central values and the lower 
bound is the value of d 2 - Since /x is rather insensitive to effects of the hnite volume, the NLO 
contact interaction with vertex —iD 2 fi^ barely differs in a hnite box. The renormalization to 
Ex = 0.2 MeV at physical pion mass then explains the similarity of the outcome for d 2 = 0 
and d 2 = —1. The renormalization condition further explains why there is no error band for 
the binding energy in the inhnite volume for mjr = 135 MeV. The contribution of the NLO 
contact interaction with coupling constant C 2 on the other hand is proportional to (7^)^ 
and since hnite volume corrections to 7 are signihcantly greater than those to /x, the error 
band for physical pion mass is predominantly determined by tq. 

The binding energy is, as expected, increasing for decreasing box size and approaching 
the inhnite volume value for large volumes. However, even at L = 20 fm hnite volume con¬ 
tributions are still above 50%. The X is signihcantly deeper bound for small box lengths 
and hnite volume corrections yield the dominating contribution to the binding energy. Be¬ 
sides the demand that pions can be included perturbatively, it is required that the binding 
momentum does not exceed the scales integrated out, which are at the order of the pion 
mass. For a volume with L > 5 fm. Ex ^ 4 MeV corresponding to binding momentum 
< 90 MeV, we expect that XEFT properly describes the dynamics of the X. 

The chiral extrapolations for hxed box size of L = 10 fm, L = 15 fm and L = 00 are 
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FIG. 7: Volume dependence of the binding energy for physical quark masses. We use the same 
notation as in Fig. 6. Note that the effective range expansion breaks down for pion masses below 
142 MeV and hence effective range results are not included. Furthermore, we renormalized the 
binding energy in the infinite volume to 0.2 MeV and thus the infinite volume result is represented 
by a single line. 





FIG. 8: Quark mass dependence of the binding energy for various box lengths. From top to bottom, 
box lengths of L = 10 fm, L = 15 fm and L = oo are accounted for. We use the same notation as 
in Fig. 6. Again, no effective range results are included since these are not valid over the whole 
range of the quark masses. 
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shown in Fig. 8. The inhnite volume results are again depicted by solid lines. The NLO 
parameters for the bounds coincide with the ones for Fig. 6. As in the inhnite volume, the 
binding energy in a hnite box shows only a moderate sensitivity to the light quark masses. 
The central values of the hnite volume belong to (i2 = 0 and tq = 0 and approach the lower 
bound for physical pion mass for reasons explained above. 


VI. CONCLUSION AND OUTLOOK 

In this work, we examined the X(3872) in a hnite volume using XEFT to NLO. We 
combined our results with chiral extrapolations, which we derived in an earlier paper [30]. 
A feature of XEFT is that NLO interactions can be included perturbatively as long as 
the expansion parameter for the inclusion of pions is sufficiently small. Based on rather 
conservative assumptions, we estimated domains for the light quark masses and, in the 
hnite volume, for the box size, where XEFT is expected to remain valid. 

In these domains, we gave explicit expressions for hnite volume corrections to the binding 
energy. On the one hand, we utilized a method, used in the inhnite volume as well, which 
can be applied for all considered values of the quark masses, even those for which the D*^ can 
decay. On the other hand, for stable we further employed an approach implementing 
the ehective range expansion, which served as a consistency check. 

Moreover, we showed that the hnite volume shift to the binding energy is fully determined 
by inhnite volume parameters and no additional input is needed. By implication this means 
that the two undetermined parameters of XEFT to NLO, denoted by tq and d 2 -, can be 
determined from lattice calculations. Here we estimated natural ranges and varied them 
within these to determine the error bands. Although there are certain values of tq and 
(i2, where the results of the two strategies mentioned above deviate and the error bars are 
possibly underestimated, over most of the natural ranges for the NLO parameters both 
methods yield consistent results. 

For all examined values of box sizes and light quark masses, we found that hnite volume 
corrections play a crucial role and yield shifts being at least in the order of the physical 
binding energy. Furthermore, over the whole natural ranges of the NLO parameters rg and 
d 2 , the D^D*^ system is bound. From these hndings, we conclude that the X should be 
observable on the lattice and already at box lengths L ~ 20 fm is expected to be a factor 
two more deeply bound than experimentally measured. 

Our analysis could be used in order to extrapolate results of lattice simulations to physical 
quark masses and inhnite volumes. In the hrst full lattice QCD study of the X(3872) in 
2013, Prelovsek and Leskovec [24] identihed a state (11 ±7) MeV below the D^D*^ threshold 
with the X for squared pion masses about four times the physical value and a spatial box 
size of 2 fm. This pion mass is clearly beyond the range of applicability of XEFT. Ignoring 
this problem and extrapolating the physical binding energy to this volume and pion mass 
using our results, we hnd the value = 270 MeV, L = 2 fm) 20 MeV which is 

roughly consistent with the lattice result. A newer lattice study from 2015 [26] arrives at a 
similar result as Prelovsek and Leskovec [24] as does the preliminary outcome from Lee et 
ah [27]. 

Among the outstanding challenges are the systematic incorporation of discretization ef¬ 
fects, an analysis of the impact of unphysical charm quark masses and an improved under¬ 
standing of the effect of cc operators [56]. Finally, an analysis of coupled channel effects, in 
particular the determination of the ^Si — ^Di mixing angle, e.g. by following the strategy 
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described for nucleons in [57], remains for future work. An analysis of the finite volume 
corrections in the Galilean-invariant version of XEFT [22] would also be interesting. 
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Appendix A: Next-to-leading order scattering diagrams 


The NLO off-shell scattering diagrams depicted in Fig. 4 read 
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where £ and £! are off-shell momenta while r] = a/— depends on the 
energy. Moreover, R = + log (tt)) is a renormalization constant, which we absorb in 

the NLO coupling constant D 2 (cf. Eq. (12b)). 


Appendix B: Next-to-next-to-leading order cancellations 

Starting from Eqs. (16) and (17), we first compute the shift of the residue up to NLO by 
equating the terms linear in the LO amplitude and obtain Eq. (18a). Plugging the shifted 
residue in Eq. (17) yields for the shift of the binding energy 

ArNlo ^ (B.l) 

1 + Si 

The dimensionless coefficient Si derives from the NLO amplitudes and is thus expected to 
be much smaller than 1. Hence, we expand the overall factor 1/(1 -|- si) in Eq. (B.l) as a 
geometrical series and obtain 

^^NLO ^ ^ 0 {s\s2)) . (B.2) 
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Taking into account that S 2 is of NLO, too, we anticipate that terms proportional to S 1 S 2 
are actually of NNLO. In fact, at NNLO there are six diagrams proportional to the LO 
amplitude squared, depicted in Fig. 9, leading to an NNLO shift for the coefficient S 2 , 
which exactly cancels out the second term in Eq. (B.2). This can be seen utilizing that 
amplitudes separate at resummed LO vertices. It is because of this cancellation that we 
omit the factor 1/(1 + si) in Eq. (18b). 

NDc\ 




FIG. 9: Scattering amplitudes shifting the coefficient S 2 in Eq. (13) at NNLO and causing a 
cancellation in Eq. (18b). 


Appendix C: Calculation of the one-pion exchange amplitude 
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The OPE amplitude ^ 0 ^^^ is depicted in Fig. 4. We begin with the unregularized expres¬ 
sion in the infinite volume 
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To transition into the finite volume we replace the spatial integration by sums over the 
allowed lattice momenta 
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At the same time we keep the contour integration over the time component since lattice 
simulations are usually performed with signihcantly larger time than spatial interval. We 
acquire 
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As a next step, we evaluate at an energy E = neglect terms proportional to 

mT,/mD and respectively, and use a tensor decomposition to replace 
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We get for the scalar amplitude 
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To renormalize, we introduce a cut-off A„ for the sum, add and subtract the inhnite volume 
loop integrals evaluated at zero energy, one regularized using a cut-off and the other one 
using PDS and take the limit —)■ cx) 
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where S{x) is dehned in Eq. (23) and S(iy){x,y) is given as 
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For the cut-off regularized integral, we expanded in /i/A 
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and found for the coefficients 


4 ^^ = 27^^ 4^^ = 27r^ (-1 + in ), 

jW = 0, = -8i7^^ 

4^^ = —27r^, = —2n‘^ {1 + in + log(4)). 
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